Shear diversity prevents collective synchronization 
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Large ensembles of heterogeneous oscillators often exhibit collective synchronization as a result of mutual 
interactions. If the oscillators have distributed natural frequencies and common shear (or nonisochronicity), 
the transition from incoherence to collective synchronization is known to occur at large enough values of the 
coupling strength. However, here we demonstrate that shear diversity cannot be counterbalanced by diffusive 
coupling leading to synchronization. We present the first analytical results for the Kuramoto model with dis- 
tributed shear, and show that the onset of collective synchronization is impossible if the width of the shear 
distribution exceeds a precise threshold. 
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PACS numbers: 05.45.Xt 

Collective synchronization is a form of self-organization in 
time that results from the interactions among a l^e hetero- 
geneous population of self-sustained oscillators [1-5]. This 
phenomenon is observed in a large variety of systems that 
range from biology to chemistry, physics and engineering (see 
e.g. [6]). For the sake of mathematical simplicity, most theo- 
retical advances in this field consider oscillators with different 
natural frequencies. Nevertheless, it is of great importance 
to know how other sources of heterogeneity influence collec- 
tive synchronization. This question has been addressed in the 
cases of heterogeneous patterns of connectivity ||7|] and inter- 
action delays [8]. 

The so-called shear (or nonisochronicity) is a crucial non- 
linear ingredient for the formation of patterns in oscillatory 
extended media l|2l|9|], as well as for the onset of complex be- 
haviors in ensembles of identical limit-cycle oscillators llioll . 
Although shear is a generic feature of oscillators, studies con- 
sidering shear diversity in ensembles of oscillators are very 
scarce and focus on a regime far from synchronization 11 1111 . 
Here we will show that distributed shear plays a key role in 
collective synchronization and may even prevent its onset. 

Mathematical formulations of collective synchronization 
usually consider as elementary oscillatory unit a normal form 
describing a system near the onset of oscillations (via the Hopf 
bifurcation), the so-called Stuart-Landau (SL) oscillator [2]: 



g = g{l- g^), 



+ 9(1-^' 



(1) 



Here the natural frequency lu determines the frequency of ro- 
tation on the attractor of radius g{t) = 1. Additionally q 
quantifies the shear of the flow, i.e. how much perturbations 
off the limit cycle modify the angular frequency 9. It is then 
usual to consider an ensemble of A^ ^ 1 globally coupled SL 
oscillators ([T]l, a mean-field version of the complex Ginzburg- 
Landau equation with disorder 111211 : 

Zj = Zj [l + i{ujj + Qj) — (1 + iqj)\zj\'^'\+K {l+ici){Z— Zj) , 

(2) 
where z^ = g^e'^^ j = 1, . . . , iV, and Z = A^^^ Y.k=i ^k- 

Previous works studying model ^ have adopted the sim- 
plifying assumption that diversity is only present in the natural 



frequencies ujj, and the shear is either absent qj — lll2Lll3ll . 
or constant qj = q il4\ . However, in a heterogeneous ensem- 
ble, either the inherent (say genetic variability for living or- 
ganisms, or tolerances for electronic circuits) or the imposed 
(e.g. experiments using coupled chemical reactors OlSll ) disor- 
der will generically be reflected in both natural frequency and 
shear terms. 

The aim of this Letter is to analyze the genuine problem 
of collective synchronization in a large ensemble of SL os- 
cillators (O with u!j and qj distributed. In our mathematical 
analysis we will assume that the oscillators are weakly cou- 
pled, i.e. \K\ is small. In such case, the dynamics of Eq. (O 
is well described by the phases only |2], 

6 J = LUj + Kqj + A:i?[sin(* ~ Sj) - q^ cos(* - 9^)]. (3) 

Here, the complex order parameter r — i?e** ~ 
N~^ J2k=i ^**'° ^^ ^ mean field and measures the degree of 
synchronization in the population. For the sake of simplic- 
ity, we assume a purely dissipative coupling, ci — 0. In 
Eq. (|3]l, the well-known Kuramoto model is recovered in the 
fully isochronous case, qj — 12l|5lll2|,[l6|], whereas the non- 
isochronous case without disorder, qj = q , corresponds to the 



so-called Sakaguchi-Kuramo to model II17I1 BlSll . 

To analyze model (jSj we adopt its thermodynamic limit 
N —> oo. Thus we drop the indices and introduce the prob- 
ability density for the phases f{9,u},q,t) il9ll . Then, the 
quantity f{9, lu, q, t) d9 du) dq represents the ratio of oscilla- 
tors with phases between 9 and 9 + d9, natural frequencies 
between lo and lo + duj, and shear between q and q + dq. The 
density / obeys the continuity equation 

dtf = -de [Uj + Kq + I [re-'\l - zq) - c.c] 1 / 

(4) 
where c.c. stands for complex conjugate of the preceding term, 
and the complex order parameter is 

/ / e'''f{9,uj,q,t)d9dujdq, (5) 

If the phases are uniformly distributed r vanishes. This state 
is customarily referred to as incoherence. Since f{9, to, q, t) 



is real and 27r-periodic in the 9 variable, it admits the Fourier 
expansion 



f{9,u,q,t) = 



2tt 



Y^ fi{u,q,t)e'' 



(6) 



/— — oo 



where // — /* ;, /o — 1, and p{uj, q) is the joint probability 
density function (pdf) of lu and q. The first Fourier mode is 
important because it determines the order parameter (|5]i: 



r*{t) 



Pi^,q)fii^,q,t) dudq. 



(7) 



Inserting the Fourier series (|6]l into the continuity Eq. dU we 
obtain the following set of integro-differential equations for 
the Fourier modes 

Kl 
dtfi - -^il{uj+Kq)fi + — [r*{l + iq)fi^i - r(l - iq)fi+i] 

18) 



Recently Ott and Antonsen (OA) found that the ansatz i20l,l21ll 

fl{u,q,t)^a{uj,q,t)', (9) 

is a particular solution of the Kuramoto model and related sys- 
tems with distributed natural frequencies. Here we also resort 
to (|9]l as it turns out to be a solution in our case if a obeys 

dta = -i{uj+Kq)a+— [r*(l + iq) - r {I - iq) a^] (10) 

The idea behind Ott and Antonsen's approach is to simplify 
an infinite set of equations — like Eq. ( fTOl l — using distributions 
that can be inserted in Eq. (|7]i and integrated via Cauchy's 
residue theorem (see below). 

In this Letter we assume that lo and q are independent ran- 
dom variables, p(w, q) = g{uj)h{q). Moreover we restrict the 
analysis to symmetric unimodal pdfs g{uj) and h{q) centered 
at cjq and q^, respectively. We start choosing g and h to be 
Lorentzian (Cauchy) pdfs. 



3(w) 



S/tt 



{uj - a;o)2 + 52 



h{q) 



7/1" 



iq-qo^+i^ 



(11) 



The integrals in (jT) can be evaluated by means of the residue 
theorem with the contour closings at infinity in the lower or 
the upper half plane of C, granted a — fi can be continued 
from real lo and q into complex lu — tOr + i^i and q = qr + iqi- 

Regarding variable a;, analyticity of a holds in the lower 
half complex oj-plane (see 1,20,1 ). As g{uj) = (27ri)^^[(w — 
ujo — i6)^^ — (w — Wo + i6)~^] has only one simple pole 
ujP = luq — iS inside this integration contour, only the value of 
a = f I at LU = Lu'P counts in the integral over lu in Eq. (|7]l. 

The integration over q in Eq. (|7]l is more intricate. We have 
to choose an integration contour such that, if a is analytic and 
|a| < 1 everywhere inside the contour at i = 0, this will hold 
for all i > 0. As a is a solution of Eq. dTol l. the analyticity of a 
at i = is preserved as t grows if a remains finite |22]. More- 
over, if a is analytic, the Cauchy-Riemann conditions imply 



dq^ |a| + i9g. |a| > 0, and the maximum of |q;| is necessarily 
located on the boundary (namely, on the integration contour). 
First of all, setting a = |a|e~*^ in Eq. dTOl l. we obtain that, 
on the real q-axis, |a| is governed by 



dt\a\ = -S\a\ + ^ Re [r*e''^(l + iq)] (l - 



(12) 



The fact that dt\a\ = —S < Oatjal = 1 guarantees that, 
if an initial condition satisfies \a{u!,q,t = 0)| < 1, this will 
hold for all t > 0. Consequently the series in Eq. (|6|l remains 
convergent. Regarding the line integral along the semicircular 
path q — |g|e''' with \q\ — > 00, Eq. ( fTOl l yields 

^ = |«|sin^-|[sin(^ + x) + |a|'sin(^-x)] (13) 
K\q\ 2 

where x{^,Q,t) = -0(^,9, t) - *(i). At \a\ = 1 Eq. (O 
becomes 



dt\a\ == (1 - i? cosx)K\q\ sini? 



(14) 



In this equation the desired relation, dt\ct\ < 0, is fulfilled in 
the lower half complex g-plane -d G (— tt, 0) only if ii' > 0, 
and in the upper half -plane -d E (0, tt) for K < 0. 

The integral over q in (|7]i can be now conveniently evalu- 
ated, and yields 



r*{t)^a{LU^LuP,q^qP,t)^a{t) 



(15) 



with qP — qa — ij for K > 0, and qP ^ qo + ij for K < 0. 
Thus, among the infinite set of equations, Eq. ( fTol i. only the 

one at (w, g) = (w^, q^) is needed: 

a = -iujPa + —{l~iqP){l-\a\^)a. (16) 

The dynamics of the radial component |a| = i? obeys 



R 



-5+-{lTl)il~R^) 



R 



(17) 



where the "^" stands for "— " if K > 0, and "+" if K < 0. 
The incoherent state, i? = 0, is always stable except above the 
line 



Kc = 



25 



if 7 < 1 . 



(18) 



At Kc a stable nontrivial solution, corresponding to a partially 
synchronized state, appears with 



R' 



K~Kc 
K 



(19) 



Equations ( fT9] l and ( fTsT i are depicted in Fig. 1(a) and com- 
pared with numerical simulations of Eq. (|2]l, i.e. using SL 
oscillators [23J. In the Kuramoto model — recovered when 
7 = — a large enough coupling strength K always results 
in partial synchronization of the population above Kc, for any 
width 5 of the frequency distribution g(w). However, here we 
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FIG. 1. Phase diagram for (a) Lorentzian pdfs giyi) and /i(g), 
Eq. dllt . and (b) Gaussian pdfs. (a) Solid line: synchronization criti- 
cal coupling, Eq. U8t . Inset: Order parameter as a function of K 
for 7 = 0.6; the solid line obeys Eq. l |19t . (b) Solid line: crit- 
ical coupling given by Eq. l l21t . In both panels, (a) and (b), the 
symbols correspond to numerical results obtained using an ensem- 
ble of SL oscillators, Eq. l[2j, with parameters: (a) TV = 40000, 
(5 = 0.01, a;o = go = |; (b) iV = 22500, a = 0.02, a;o = go = 0. 
Parameters cjj and qj were deterministically selected to represent 
p(a;, q) and averages were done over time. 



find that the width 7 of the shear distribution h{ci) has a more 
severe effect on the synchronization transition. If 7 > 7d = 1, 
synchrony disappears for all K and (5, and incoherence be- 
comes the only stable state. It is noteworthy that this is a 
collective phenomenon caused by the presence of distributed 
shear, and it has no counterpart in the case of two coupled SL 
oscillators |241 |ll,i25ll . 

To confirm the generaUty of these findings for other pdfs 
112611 . next we follow f 19] and perform the linear stability anal- 
ysis of Eqs. dHJ about the incoherent state /j^o — 0. We find 
that the only potentially unstable modes are Z = ±1. Inserting 
/i(a;, q, i) = 6(0;, q)e^' into Eq. (O, the discrete spectrum of 
eigenvalues A can be obtained by virtue of a self-consistency 
argument. This yields the integral equation 



2 



l + iq 



X + i{uj + qK) 



g{uj)h{q) cLu dq. 



(20) 



The border of unstable incoherence Kc, is found imposing 
Re(A) -^ 0+. If g{uj) and h{q) are Gaussian functions with 
variances a^ and i'^, respectively, and h{q) has zero mean 
(go = 0), the critical coupling can be explicitly obtained: 



Kr 



-TT + 8z^2 _|_ ^7r(7r -I- 16t/2) 



(21) 



This function is plotted in Fig. 1(b) and compared with the re- 
sults of numerical simulations of Eq. (|2]). Remarkably, we find 
again a threshold for the dispersion of h{q), Vd = •\/7r/2 — 
1.253 . . ., above which incoherence is stable for all K. 

We have found that the divergence of Kc does not only 
exist for Lorentzian and Gaussian distributions, but for any 
symmetric unimodal distribution. This divergence occurs at a 
shear diversity that is conveniently expressed in terms of the 



peak value h{qo). If qo = 0, the divergence occurs at 

1 



hiO) = -. 

TT 



(22) 



Otherwise, if h{q) is not centered at zero, Kc also diverges 
at a certain value of h{qQ) — hd, but a simple distribution- 
independent formula like ( |22] | does not exist 12711 . 

Our numerical calculations of Eq. (|2]) using Gaussian pdfs 
also reveal that incoherence and synchronization coexist in a 
region with large K/a [light-shaded region, Fig. 1(b)]. Note 
that this region is not present in the Lorentzian case [Fig. 1(a)], 
in spite of the similar bell-shaped form of these two distribu- 
tions. Figure 1(b) suggests that the destabilization of inco- 
herence may also occur through a subcritical bifurcation for 
certain pdfs. To elucidate the supercritical or subcritical char- 
acter of the synchronization transition, we carry out a self- 
consistency analysis "a la Kuramoto" ||2|, [l2|, [l6|, [iTD in the 
limit of large coupling and/or very small frequency disper- 
sion, i.e. Kg{ujQ) ^ 1. After going into a rotating framework 
Ldot, rescaling time t — ^ K^^t, and neglecting the 



ojj/K term, we approximate Eq. (O by: 

dj = q.j + i?[sin(* - Oj) - q.j cos(* - 6*-,)]. 



(23) 



Hereafter we assume q^ = 0. In a partially synchronized state 
the population splits into two groups, the synchronized (or 
locked) subpopulation with \q\ < q^ax = i?/\/l — R^, and 
the desynchronized (or drifting) one with \q\ > g„iax- Both of 
subpopulations do contribute to the order parameter: 



R 



' ds 



(24) 



where we have chosen a reference frame where ^ = and 
the brackets denote averages over each subpopulation. We can 
now make an expansion in powers of R for each contribution. 
Up to cubic order we obtain: 
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R'^-h{0) + R^ 



MO) 



h"{0) 



R'-h{0) 



i?3 r 
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h'{q) 



dq 



where we have assumed h{q) is twice differentiable at the ori- 
gin. Inserting these expansions into Eq. (l24l l. and discarding 
the trivial solution i? = we find that R follows asymptoti- 
cally a square-root dependence of the form: 



R 



1 - Trh{0) 



J 



where 



-^ 



7:h"{0) 



h'iq) 



dq 



(25) 



(26) 



is evaluated at h{0) = tt^^. The sign of J determines the ori- 
entation of the bifurcating branch in Eq. ( |25] |, as showed in the 
conjectured bifurcation scenarios in Fig. 2(a,b). When J > 0, 
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FIG. 2. (a,b) Sketch of the possible bifurcation scenarios for large 
K values, depending on the sign of J, Eq. l l26b . Solid (dashed) 
lines indicate stable (unstable) solutions, (c) Bifurcation diagram for 
Lorentzian pdfs at K/5 = 40, obtained from Eqs. l |19t and l llSt . 



as for the Gaussian distribution (J = 



If = 0.261...), 



a partially synchronized solution branches off from incoher- 
ence subcritically. This is in concordance with the numerical 
results in Fig. 1(b). The scenario for J < 0, Fig. 2(b), is also 
followed by pdfs with a non-differentiable maximum of type 
h{q) = h{0) + h'{Q~^)\q\ + • • • , like the triangular or Laplace 
distributions (formally J = — cx)). It is important to note that 
bistability incoherence-synchronization is found irrespective 
of the sign of J, because Eq. (|23T l has a stable fixed point at 
6'j = ^ (i? = 1) that persists in the form of a solution with 
R near 1, provided g{Lu) has a small dispersion (or K is large 
enough). 

The case of Lorentzian h{q) is quite peculiar. On the one 
hand, J vanishes for this pdf, what is consistent with the in- 
finitely abrupt transition predicted by the OA ansatz in the 
limit K — > oo, see Fig. 2(c). On the other hand, according to 
our numerical simulations, the synchronized solution R = 1 
of Eq. ( |23] | does not persist for 7 > 1 if g{uj) is not a delta 
function. This indicates that, for heavy-tailed h{q), the term 
ujj / K neglected in Eq. (|23] | may become relevant and destroy 
the synchronized solution beyond a certain critical value of 
h{,qo)-\ 

In summary, this Letter uncovers the effect of shear diver- 
sity on the collective synchronization of globally coupled os- 
cillators. We have obtained the first analytical results for the 
Kuramoto model with distributed shear (|3]l. They show that, if 
shear is widely distributed, incoherence is always stable and 
for some distributions — such as the Lorentzian one — syn- 
chronization is impossible. The techniques used here can be 
readily applied to a number of extensions of the model ^, 
such as considering other distributions p{uj^ q), periodic or 
stochastic driving, time delays, or networks and communities 
of oscillator populations. 
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